Bounded dynamics of finite PT-symmetric magnetoinductive arrays 
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We examine the conditions for the existence of bounded dynamical phases for finite PT-symmetric 
arrays of split-ring resonators. The dimer (N = 2), trimer (N = 3) and pentamer (N = 5) cases 
are solved in closed form, while for N > 5 results were computed numerically for several gain/loss 
spatial distributions. It is found that the parameter stability window decreases monotonically with 
the size of the array. 



The study of PT-symmetric systems has attracted a 
lot of attention during the past few years. In these sys- 
tems, the effects of loss and gain can balance each other 
and, as a result, give rise to a bounded dynamics. These 
studies are based on the seminal work of Bender and 
coworkers [H [1| who showed that non-hermitian Hamil- 
tonians arc capable of displaying a purely real eigen- 
value spectrum provided the system is symmetric with re- 
spect to the combined operations of parity (P) and time- 
reversal (T) symmetry. For one-dimensional systems the 
PT requirement leads to the condition that the imaginary 
part of the potential term in the Hamiltonian be an odd 
function, while the real part be even. The system thus 
described can experience a spontaneous symmetry break- 
ing from a PT symmetric phase (all eigenvalues real) 
to a broken phase (at least two complex eigenvalues), 
as the gain/loss parameter is varied. To date, numer- 
ous PT-symmetric systems have been explored in several 
fields, from optics [S3, electronic circuits Q, solid-state 
and atomic physics [3.1 lOf. to magnetic metamaterials 11} . 
among others. The PT symmetry-breaking phenomenon 
has been observed in several experiments (? . Tl2 . 13j ]. 



Magnetic metamaterials, on the other hand, consist of 
artificial structures whose magnetic response can be tai- 
lored to a certain extent. A common realization of such 
system consist of an array of metallic split-ring resonators 
(SRRs) coupled inductively l44l6|| . This type of system 
can feature negative magnetic response in some frequency 
window, making them attractive for use as a constituent 
in negative refraction index materials [13] ■ A common 
problem with SRRs is the heavy ohmmic and radiative 
losses. One possible solution is to endow the SRRs with 
external gain, such as tunnel (Esaki) diodes 18, [l9| to 
compensate for such losses. 

In the case of some one-dimensional coupled discrete 
systems, such as a harmonic oscillator array, it has been 
observed that in the limit of an infinite size array, the sys- 
tem belongs to the broken PT phase, i.e., there are com- 
plex eigenvalues making the dynamics unbounded [HI [20| . 
In this work we examine the case of short SRRs arrays, 
and determine the parameter window inside which the 
system exhibits a bounded dynamics, and how this win- 



dow decreases with the size of the system. 

Let us consider a simple model of a magnetic meta- 
material consisting of a finite one-dimensional array of 
split-ring resonators including gain/loss terms: 
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where A is the magnetic interaction coefficient or cou- 
pling, and 7„ positive (negative) denotes a ring with 
loss (gain). In order to satisfy the requirements for PT- 
symmetry, the spatial distributions of the gain/loss must 
be odd, 7_n = — 7«. In this work we will focus in binary- 
like systems with two gain/loss terms and thus examine 
arrays of the form . . . - 71, -72, -71, 0, 71, 72 , 71 , . . ., for 
arrays with an odd number of rings. For arrays with 
even number of rings the distribution of gain/loss is of 
the form ... - 71, -72, -71, 71, 72, 7i, • • •■ This distinc- 
tion is only meaningful for small arrays and disappears 
for system of infinite size. Hereafter, and without loss of 
generality, we will focus on arrays with an odd number 
of sites (except for the dimer case). Results for the case 
with an even number of sites are similar. Since the values 
of 71 and 72 are arbitrary, the gain/loss distribution thus 
introduced allows for many interesting cases to be exam- 
ined. In particular we will focus on three cases. The first 
one is 71 = 7, 72 = —7, giving rise to the distribution 
. . . 7, — 7, 7, — 7, 0, 7, — 7, 7, — 7. . .. The second one is 
71 = 7 and 72 = 0, which gives rise to the distri- 
bution ... 0, —7, 0, —7, 0, 7, 0, 7, 0, . . .. Another interest- 
ing case is the one with 71 = 7 = 72 that gives 
. . .—7, —7, —7, —7, 0, 7, 7, 7, 7. . .. For this last case, we'll 
see that in spite of the concentration of loss and gain on 
opposite sides, the dynamics does possess a stability win- 
dow for finite arrays lengths. 

Dimer case (N = 2): In this case, the only possible case 
is 7, —7. The dynamical equations read: 
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T ^(qi+Xq 2 ) +7-r<7i+<Zi = (2) 
dr z dr 

d 2 d 

— (q 2 + Xqi) - 1— qi + q 2 = (3) 
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Figure 1: Dimer array. Left: Stability region (shaded) in 
gain/loss-coupling space. Right: Mode frequency squared as 
a function of the gain/loss parameter for A = 0.3. 
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Figure 2: Dimer array. Mode phase as a function of the 
gain/loss parameter. Solid: Q ++ mode. Dotted: Q + ~ mode. 
Short dash: mode. Long dash: f2 mode. (A = 0.33) 



We look for stationary modes qi i2 (r) = qi t2 cxjp(iQr). 
This leads to the equations, 
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The condition of the vanishing of the determinant of this 
linear system leads to a quadratic equation for SI 2 , with 
solutions: 
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We denote the four solutions as S1 ++ ,S1 H , SI h and 
S7 . The stable phase (unbroken PT symmetry) cor- 
responds to the cases where SI is a real quantity. From 
straightforward examination of Eq.(6), one concludes 
that the stability window in 7- A space is given by the area 



under the curve 7c = y 2(1 — vl — A 2 ), for < A < 1. 
Outside 7 = 7c (A), the system is unstable. Figure Q] 
shows the stability region and also the square frequencies 
as a function of the gain/loss parameter. Due to symme- 
try considerations, only the positive 7- A sector needs to 
be explored. 

From Eq.(4) ( or Eq.(5)) and Eq.(6), it is easy to obtain 
I92/91I = 1 & n d thus, qi and q 2 differ by a phase only. 
We have four branches for the phase, corresponding to 
each of the four solutions. Figure [2] shows the phase of 
all solutions as a function of the gain/loss parameter. 



Trimer case (N = 3). 
has the form —7,0,7. 
have the form 



Here the gain/loss distribution 
The stationary state equations 
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with eigenvalues 



O 
O 



±1 



± 



2 - 7 2 ± y/-4 7 2 + 7 4 + 8A 2 
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The condition that all Q be re al leads to th e condi- 
tions A < l/y/2 and 7 < j c = 8A 2 . Figure 
[3] shows the window of real eigenvalues in gain/loss- 
coupling space. It is qualitatively similar to the dimer 
case, but the allowed coupling interval is smaller. The 
figure also shows the square of the real frequency as a 
function of gain/loss, for a given coupling value. 

Pentamer case (N = 5): Here the gain/loss distribu- 
tion can have three possible forms: 7,-7,0,7,-7, or 
0,-7,0,7,0, or — 7, — 7, 0, 7, +7. We will focus on the 
first case since is more amenable to an exact form so- 
lution (Numerical results show that the other two cases 
display similar behavior). The stationary state equations 
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Figure 3: Trimer array. Left: Stability region (shaded) in 
gain/loss-coupling space. Right: Mode frequency squared as 
a function of the gain/loss parameter, for A = 0.3. 
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Figure 4: Pentamer array. Left: Stability region (shaded) in 
gain/loss-coupling space. Right: Mode frequency squared as 
a function of the gain/loss parameter, for A = 0.3. 



have the form 
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with eigenvalues 
fi = ±1 
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The condition that all tl be r eal leads to the conditions 
A < 1/V3 and 7 < 7c = ^2(1 - y/T- A 2 ). Figure H 
shows the stability window in gain/loss-coupling space, 
as well as the square frequency as a function of gain/loss, 
for a fixed coupling value. As we can see, the stabil- 
ity window is substantially smaller than the one for the 
dimer and trimer cases. 

Short chains (N > 5): Let us consider now the case 
of finite arrays, where the stationary-state equations are 
given by 

- n 2 [q n + X(q n+1 + g n _x)] + (1 + in ln )q n = 0.. (20) 

The {7n} distribution we consider has the general form: 

. . . — 71, —72, —71, 0, 71,72,71, We will compute the 

relevant eigenvalues numerically from the vanishing of 
the determinant of linear system ([20]) and will focus on 
the parameter values in gain/loss - coupling space where 
the eigenvalues are purely real, thus denoting a bounded 
dynamical regime. 




Figure 5: Stability regions (area under curve) in gain/loss- 
coupling space for several array lengths, for the case 
. . . , -7, +7, -7, , +7, -7, +7, . . . . 



Case a. We start with the case 71 = 
—72 = 7 that gives rise to the distribution 

. . . , —7, +7, —7, O, +7, —7, +7, Some results 

are shown in FigfS] It is clear that, as the size of 
the array increases, the stability region shrinks, and 
disappears altogether around N = 20. This is consistent 
with a previous result 
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stating that in the infinity 
size limit, the dynamics is always unstable (that is, the 
system belongs to the broken PT phase). 
Case b. Now we take 72 — > and 71 = 7, that is, the 

distribution . . . , 0, —7, 0, —7, O, +7, 0, +7, 0, +7, 0, 

Notice how the gain and loss portions are now separated 
and on each side they are rather diluted. The stability 
phase diagrams for this case (not shown) are qualita- 
tively similar to the previous case, although, for a given 
array size, the stability windows are smaller. 
Case c. Now we take 72 = 71 = 7. The gain/loss distri- 
bution is 

. . . , —7, —7, —7, —7, O, +7, +7, +7, +7, . . . Now the gain 
and loss on each side are densely populated and the re- 
sulting area of the stability windows (not shown) while 
qualitatively similar to the previous cases, drop even 
faster. 

Figure [6] shows a summary of the results obtained for 
the size of the stability window as a function of the array 
length, for the three cases considered. We clearly see that 
the stability region decreases rather abruptly with N, 
and that for a given N, we have area a < area b < area c, 
for N > 5. 
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Figure 6: Normalized size of area in gain/loss-coupling space 
with purely real eigenvalues as a function of array size for the 
gain/loss distributions "a" (squares), "b" (circles), and "c" 
(triangles). 



Conclusions. We have examined the dynamics of finite, 
one-dimensional PT-symmctric arrays of split-ring res- 
onators, that constitute the simplest model of a mag- 
netic mctamaterial, and have focusscd on the conditions 
in parameter space where this dynamics is bounded, i.e., 
all eigenvalues are real. The dimcr, trimer and pentamcr 
cases were solved in closed form, while for larger but finite 
arrays, results were obtained numerically. It was found 
that in general, the stability region decreases abruptly 
with an increase in system size. This tendency seems 
generic for gain/loss distributions based on two gain/loss 
parameters, and is in agreement with recent results for 
infinite discrete arrays [ill l20l]. 
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